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Abstract: We study the origin of subleading soft and collinear poles of form factors 
and amplitudes in dimensionally-regulated massless gauge theories. In the case of form 
factors of fundamental fields, these poles originate from a single function of the coupling, 
denoted G{as), depending on both the spin and gauge quantum numbers of the field. We 
relate G{as) to gauge-theory matrix elements involving the gluon field strength. We then 
show that G{as) is the sum of three terms: a universal eikonal anomalous dimension, a 
universal non-eikonal contribution, given by the coefficient Bs{as) of 5{l — z) in the collinear 
evolution kernel, and a process-dependent short-distance coefficient function, which does 
not contribute to infrared poles. Using general results on the factorization of soft and 
collinear singularities in fixed-angle massless gauge theory amplitudes, we conclude that 
all such singularities are captured by the eikonal approximation, supplemented only by 
the knowledge of Bs{as). We explore the consequences of our results for conformal gauge 
theories, where in particular we find a simple exact relation between the form factor and 
the cusp anomalous dimension. 
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1. Introduction 

The structure of soft and collinear singularities in pcrturbative QCD has been studied in 
depth for decades, uncovering a pattern of exponentiation dictated by gauge invariance 
and factorization. The prototype amplitude for these studies is the electromagnetic form 
factor of a colored particle, whose simple color structure and renormalization properties 
make it an ideal laboratory for isolating and evaluating long-distance contributions. In 
fact, the color-singlet QCD form factors of quarks and gluons can both be expressed in an 
elegant exponentiated form, each in terms of only two functions of the running coupling. 
In this paper, we will extract and interpret the universal functions that control subleading 
soft and collinear poles for the form factors. The same functions, as we will see, also 
control collinear poles for the full class of dimensionally-regulated fixed-angle scattering 
amplitudes. 

Color-singlet parton form factors are the simplest amplitudes exhibiting the double 
logarithmic ('Sudakov') behavior characteristic of gauge theories in the massless limit. 
Following the early studies in the abelian theory [1], which were performed at leading 
logarithmic (LL) accuracy, the form factors of non-abelian gauge theory were shown to 
exponentiate to arbitrary logarithmic accuracy [2, 3, 4, 5]. Exponentiation occurs because 
the form factor obeys an evolution equation, which in turn is a consequence of factorization 
and gauge invariance [6]. Solving the evolution equation yields an especially transparent 
answer [7] if one employs dimensional regularization as an infrared regulator, as is routinely 
done in finite-order perturbative calculations. 
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Dimensional regularization, in this context, displays several remarkable features, going 
well beyond its properties of preserving gauge symmetry and simplifying calculations in 
massless theories. When performing a resummation, in fact, dimensional regularization 
expresses the solution to the appropriate evolution equation in terms of the d-dimensional 
running coupling, which vanishes in the infrared for d > 4 as a consequence of dimensional 
counting. On the one hand, this allows one to solve the equation in terms of a simple initial 
condition, since all perturbative contributions to the amplitude vanish as a power of the 
hard scale for d > 4; as a consequence, the resummed amplitude can be directly compared to 
finite-order Feynman diagram calculations. On the other hand, the d-dimensional running 
coupling in general displays a Landau pole with a nonvanishing imaginary part, which in 
turn allows an explicit evaluation of the resummed amplitude in terms of analytic functions 
of the coupling and the dimension d [8] . 

Fixed-angle scattering amplitudes also have poles in dimensional regularization, and 
more generally double-logarithmic infrared enhancements at high energy. Evolution equa- 
tions for on-shell high-energy scattering were developed first for theories without gauge 
bosons [9], while the infrared divergences associated with external colored particles re- 
mained an obstacle [10]. The recognition that soft and collinear (virtual) radiation can 
be factorized in a universal way from a hard QCD process (and specifically the result 
that soft radiation can be factorized from harder collinear radiation [11]) made it possible 
to generalize the exponentiation of the form factor to amplitudes with multiple colored 
legs [12, 13, 14]. 

Any gauge theory scattering amplitude can be treated as a vector in the space of 
available color configurations [15], 




= . (1.1) 

Here the vector associated with a scattering process with flavor structure labeled 

by [f], is represented by the coefficients A^^', in a basis defined by a set of N^"^ color tensors 
c^; /9j are particle momenta "pj rescaled by a hard scale Q, for example as 'Pj = {Q/-\/2)/3j. 
When the amplitude describes fixed-angle scattering, which we represent as 

i-Pl+P2^P3 + ■■■+Pn+2, (1-2) 

each color component can be factorized into a product of 'jet' functions, describing 

the virtual color-singlet evolution of each external hard particle due to collinear radiation, 
times a 'soft' function organizing the effects of long- wavelength radiation [12, 13, 14], 

•A^? ^,a,(/.2),e) = n (^,a,(/x2),e) ^, a,(/.2), 

xi?f (/3„^,^,a,(;.^)) , (1.3) 
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where Q' plays the role of a factorization scale separating soft and coUinear momenta. The 

[fi 

soft function, Sjj^, is a matrix in the vector space spanned by the color tensors c^. It acts 
on a vector of finite coefficient functions describing the effects of highly virtual particles. In 
this context, the form factor continues to play an important role: one can always, in fact, 
choose the factorization scheme so that the 'jet' functions for a generic amplitude may be 
identified with the square roots of the form factors of the corresponding hard partons. 



The soft matrix Sj]j in Eq. (1.3) is responsible only for purely infrared single poles. For 
this reason, it can be computed using the eikonal approximation, in which hard partonic 
lines are replaced by Wilson lines. This fact leads to the important conclusion that the 
non-eikonal, hard collinear singularities in an arbitrary QCD amplitude, for fixed-angle 

scattering, can be completely organized in terms of the color-singlet form factors of quarks 
and gluons. Our discussion below will focus on the non-leading pole structure of the form 
factors. 

Most of the results described above were originally derived with phenomenological ap- 
plications in mind: long-distance singularities in the amplitudes, in fact, are the source 
of logarithmic enhancements in infrared- and collinear-safe cross sections near kinematic 
boundaries, which can have a sizable impact on perturbative predictions and often need 
to be resummed to all orders. Furthermore, the universal structure of long-distance sin- 
gularities at fixed order [16] provides an important test for perturbative calculations, and 
is an essential ingredient in the construction of subtraction schemes, which are necessary 
to compute finite jet cross sections and event-shape distributions. Non-leading logarith- 
mic enhancements were studied in deep-inelastic scattering and vector boson production 
in [17, 18, 19, 20]. These papers uncovered the same pattern for subleading enhancements 
that we will identify for poles in fixed-angle scattering. We will return to these studies in 



Beyond immediate phenomenological consequences, our analysis has a wide range of 
applicability, although derived within the context of QCD. It is based on universal prop- 
erties of quantum field theories and of gauge theories in particular, and encodes general 
information about their long-distance behavior. In recent years, in fact, QCD results on 
the infrared structure of amplitudes have been applied to supersymmetric theories, and 
in particular to the maximally supersymmetric A/" = 4 super- Yang-Mills theory, which 
is of great theoretical interest because of its connections with string theory through the 
anti-de-Sitter-space-conformal-field-theory (AdS-CFT) correspondence [21]. 

This correspondence states that the strong-coupling, planar (large limit of = 4 
super- Yang-Mills theory admits a simple description in terms of solutions of a classical 
gravitational theory, or of strings moving in a weakly curved background. On the other 
hand, the quantum conformal invariance of the theory, which has a vanishing /? function to 

^In Refs. [13, 14] the functions r'*' were denoted as Al''^*', as in the notation for amphtudes above. 




(1.4) 



where F^ 



(Q^) represents the color-singlet form factor^ for parton i. 



Sec. 5.2. 
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all perturbative orders, implies a drastic simplification in the all-order perturbative resum- 
mation of infrared and coUinear singularities. Together, these two observations suggest that 
perturbation theory for planar A/" = 4 super- Yang- Mills amplitudes should have remarkable 
properties. 

Indeed, motivated by a surprising iterative property of the two-loop four-gluon scat- 
tering amplitude [22, 23], an exponential form for the all-loop scattering amplitude in 
planar = 4 super- Yang-Mills theory was proposed [24] . This conjecture was constructed 
to be consistent with the all-order structure of soft and collinear divergences derived in 
Refs. [7] and [13]. At the level of finite terms, it was tested successfully for the three-loop 
four-gluon amplitude [24] and two-loop five-gluon amplitude [25, 26]. For four gluons, it 
received striking confirmation from the work of Alday and Maldacena [27] , who computed 
the strong-coupling planar limit directly from the AdS-CFT correspondence, in dimen- 
sional regularization, finding the same exponential form in this limit. The conjecture is 
expected to hold for five gluons; recently, however, it has been found to break down for six 
gluons [28] . Such a breakdown was anticipated at strong coupling by an observed inconsis- 
tency in a particular kinematic limit with a very large number of gluons [29] , and at two 
loops by an analysis of the high-energy limit of the six-gluon amplitude [30]. 

Irrespective of the form of the finite terms, one can also extract from the Alday- 
Maldaccna solution the leading strong-coupling behavior of the same quantities that gov- 
ern infrared evolution in perturbation theory [27, 31]. The leading (double) poles in the 
exponent for the amplitude, for example, are governed by the cusp anomalous dimen- 
sion [32, 33, 34, 35]. The AdS-CFT correspondence can then be tested at the level of 
elementary fields, by looking for consistency between the perturbative series at weak cou- 
pling and the strong-coupling limit. 

In the case of leading poles, the comparison can be made very precise because the cusp 
anomalous dimension has a well-defined non-perturbative definition in field theory, and a 
wealth of information is available about its properties. For planar A/" = 4 super- Yang-Mills 
theory, it has been computed perturbatively up to four loops [24, 34, 36, 37, 38]. A strong- 
coupling expansion can be derived from string theory, in which the first three terms have 
now been computed [39, 40, 41]. Most remarkably, one can exploit the observed integrability 
properties of the theory to construct an integral equation [42] , whose solution extends over 
all values of the couphng, and is in precise agreement with all four weak-coupling and all 
three strong-coupling [43] coefficients. 

Subleading poles for scattering amplitudes are less well understood. In general, they 
cannot be expressed completely in terms of eikonal amplitudes, although we will see that 
eikonal amplitudes still play an important role in characterizing them. It is worth remarking 
that there has been considerable renewed interest of late in the perturbative study of eikonal 
amplitudes, or Wilson loop expectation values, inspired by their role in the strong-coupling 
approach of Alday and Maldacena. A close, but still not fully explained, relation between 
Wilson loop expectation values and the maximally-helicity-violating scattering amplitudes 
in A^ = 4 super- Yang-Mills theory has been uncovered at one loop [44] and verified at two 
loops [28, 45, 46]. 

It is the purpose of our paper to provide a more precise characterization of the 
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sublcading-pole singularities in a massless gauge theory such as QCD or = 4 super- 
Yang-Mills theory. We begin by reviewing briefly in Section 2 the known results about 
the Sudakov form factor, which we present in a simplified form and then apply to the 
case of a conformal theory. In Section 3, we revisit the standard factorization [47] in the 
context of dimensional regularization. This allows us to give explicit operator expressions 
for the functions that control the poles of the form factors, and thus the collinear poles 
of all other fixed-angle amplitudes, via Eq. (1.3) above. In Section 4, we derive evolution 
equations for these operators, and identify their anomalous dimensions. This leads to an 
explicit expression, in terms of these anomalous dimensions, for the function G{as), which 
determines the subleading poles of the form factor. In Section 5, we construct another ex- 
pression for G{as), by relating collinear singularities of the form factor, for a given parton 
species, to the virtual contributions to the corresponding parton distribution. This leads 
us to identify explicitly the only non-eikonal long-distance contribution to G{as), which is 
given by the virtual term of the diagonal Altarelli-Parisi splitting function for the chosen 
parton, a fact that was pointed out at finite perturbative order in Refs. [48, 49]. Eikonal 
contributions to G{as), on the other hand, are related to the function responsible for soft 
single logarithms in threshold resummation for the DrcU-Yan process [50]. As shown in 
Section 5 and in the Appendix, certain additional short-distance contributions to G{as) 
arc given entirely by running-coupling effects. These contributions do not give rise to poles 
in scattering amplitudes, and are proportional to e in a conformal gauge theory such as 
A/" = 4 super- Yang-Mills theory. 

We hope that our results will be helpful both in perturbative QCD studies, where a 
detailed knowledge of long-distance singularities to all orders is of direct phenomenologi- 
cal relevance, and in order to further our understanding of superymmetric gauge theories, 
where the striking discoveries of these recent years are beginning to map a precise con- 
nection between perturbation theory and strong coupling, possibly on the way to exact 
results. 

2. Gauge Theory Form Factors 

Let us begin by reviewing briefly the known results concerning the color-singlet form factors 
of massless colored particles. For a quark, one can define the form factor as a matrix element 
of the electromagnetic current. In the time-like case, for example, one can write 



Gluons do not couple directly to the electromagnetic current J^, but their form factor can 
be defined analogously as a matrix element of a gauge-invariant operator. A typical and 
useful example is the coupling of gluons to colorless scalar particles (such as the Higgs 
boson) through an effective vertex constructed by integrating out a heavy fermion loop 
(the top quark in the standard model). In this case the form factor is defined through the 
coupling to the operator — C//Tr [Gf^,i,G^'^]/2, where G^i, is the Yang-Mills field strength and 
Ch is a matching coefficient containing the dependence on the mass of the heavy fermion. 




(2.1) 
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Note that in this case the effective operator couples directly to soft gluons, but this does 
not change the structure of infrared singularities at leading power. The reasoning below 
therefore applies equally well to quarks and to gluons. 

As was shown in Refs. [2, 3, 4, 5], and reviewed in Ref. [47], the momentum depen- 
dence of the form factor is determined by a simple evolution equation. The equation is a 
consequence of the factorization of soft and coUinear modes from highly virtual exchanged 
particles and from each other [6] , which in turn arises from the loss of quantum- mechanical 
coherence for processes occurring at widely separated scales. It can be proven by making 
use of Ward identities. In dimensional regularization, with d = 4 — 2e and e < in order 
to regulate mass divergences in the renormalized theory, the evolution equation takes the 
form 



d 



log 



K (e, as{ii^)) + G ( -j^,as{i^'^), e 



(2.2) 



where the function ^(e, a^) is a pure counterterm, while the function G{(,'^ , ag, e), which 



carries the momentum dependence, is finite as e 
(RG) invariance of the form factor implies that 



0. Furthermore, renormalization-group 



d 



P{e, as 



d 



da. 



G 



, a. 



-lM^+/?(6,a.) — 



(2.3) 

which is how the cusp anomalous dimension 7ii-(Q;s) comes into play in this context. 

In order to solve Eq. (2.2), one needs to introduce the c/-dimensional running coupling, 
which solves the RG equation 

da. 



/3(e, as) = -leas + &{pL^ , 



where ${as) is the usual four-dimensional /3 function. 



(5{as 



a 
'2tt 



(2.4) 



(2.5) 



n=0 



with 6o = (IIC'a 
Eq. (2.4) is. 



2n/)/3 in our normalization. At the one-loop level, the solution to 



a 



LVMo 



1 - 



6o 



(2.6) 



which clearly reduces to the usual four-dimensional result as e — >^ 0. Since the solution is RG 
invariant, we may use for the running coupling the simplified notation a{fi^,e) whenever 
we do not need to adopt a specific boundary condition. 

At tree level (and to all orders in a conformal theory) a scales as a power of /v,, 
a{iJ?,e) ~ fi'^'^as- The running coupling thus vanishes in the infrared, as expected above 
the critical dimension d = A. Given the RG invariance of the form factor, this has the 
important consequence of providing us with a simple initial condition for Eq. (2.2), 



r(0,a,(M2),€) =r(l,a(0,e),e) = 1 



(2.7) 
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It is now straightforward to integrate Eq. (2.2), obtaining 




K{€,as{li^)) (2.8) 



where we integrated along the negative real axis to emphasize that the function G is real 
for negative Q^. 

The term proportional to K{e, as) in Eq. (2.8) has an apparent unregulated singularity 
due to the integration down to = 0: in fact, this term cancels exactly the ^^-independent 
terms arising from the integration of 7A'(a). This can be shown to all orders because the 
function K{e,as) is completely determined, through Eq. (2.3), by the coefficients of the 
pcrturbative expansions of and of the /? function [8]. As a consequence, all poles in 
Eq. (2.8) arise from integrations over the scale of the running coupling in the infrared 
region. This cancellation can be made explicit by considering the RG equation for the 
counterterm function K{e,as)-, Eq. (2.3). Since K has no explicit scale dependence, one 
can write 

d d 
^d/j,^^^' «s) = /?(e, "s) ^^(^' <^«) = -7K{as) ■ (2.9) 

Using, once again, the vanishing of the running coupling in the infrared, one has the 
boundary condition K{iJ, = 0) = 0, so that Eq. (2.9) integrates to 

K{e,as{fi')) = -^j^ ^7K(«(A^6)). (2.10) 

We now have two terms in Eq. (2.8) involving double scale integrals of the cusp anomalous 
dimension. They both diverge, but using Eq. (2.10), exchanging orders of integration, and 
choosing n = Q they can be re-expressed as 

r (Q^ e) = exp 1 1 J'"^^ ^ [g( - 1, a (e^ ^) , ^) - ^ 7x (a {f, e) ) log 

(2.11) 

In QCD, Eq. (2.11) has phenomenological as well as theoretical interest. In fact, as de- 
scribed in Ref. [13], it is one of the building blocks for the analysis of mass singularities in 

general multiparton amplitudes at finite perturbativc orders, which in turn is of relevance 
for the calculation of infrared-safe observablcs at high-energy colliders. We observe at this 
point that the function G in Eq. (2.11) not only generates next-to- leading poles in e at 
each order in a^, but it also serves as a complete infrared-safe coefficient function for the 
exponentiation of such poles, as well as finite parts, in the singlet form factor. 

In the present context, it is interesting to notice that Eq. (2.8) drastically simplifies 
in a conformally-invariant theory such as = 4 super- Yang-Mills theory. We emphasize 
again that the evolution equation, Eq. (2.2), is a consequence of gauge invariance and 
factorization, and hence holds for supersymmetric extensions of QCD, with, of course. 
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different but related functions K and G. For a conformal tlieory, /?(«<,,) = 0; tlie coupling 
then runs according to its mass dimension in d = A — 2e, so that q;s(A^)A^'^ = as{lJp)lJp''', 
as a consequence, all integrals in Eq. (2.8) can be performed trivially in this case [24]. 
Expanding the anomalous dimensions as 



oo 

E 

n=l 



(n) 

Tk 



n=l 



G(")(e) 



(2.12) 



one finds several remarkably simple results. First of all, the counterterm K(e, a^) has only 
simple poles, and is easily expressed in terms of the perturbative coefficients of jk, as 



K{e,as) 



n=l 



(n) 

Tk 
2ne 



(2.13) 



Next, one observes that the logarithm of the form factor has only double and single poles in e 
to any order in perturbation theory, in contrast to the situation in QCD, where the running 
of the coupling generates poles up to e~P''^ at order af . In fact, one finds explicitly [24], 



log 



T 

n=l 



^ oo 

-T 

n=l 



TT 



(n) 

Tk 
2n2e2 



+ 



G(")(e)' 



ne 



(n) 

Tk 
2n2e2 



+ 



G(")(e) 



ne 



(2.14) 



displaying, as expected, exact RG invariance. 

As a final remark, it is interesting to construct, in a conformal theory, the ratio of the 
time-like to the space- like form factor. This ratio was studied for QCD in Ref. [7]. In that 
case it is of phenomenological relevance, since it enters the resummed expression for the 
Drell-Yan cross section in the DIS factorization scheme [51, 52, 53]. In the conformal case 
the analytic continuation can be performed explicitly, and one finds 



log 





i 


[r(-Q2)J 


= 2^ 



K{e,as{Q^))+G{-l,as{Q^),0) +^jk {a,{Q^)) +0{e) . (2.15) 



As observed in Ref. [7] in QCD, all poles in the ratio are given by an infinite phase, which 
in this case is simply related to 7^ via Eq. (2.13); the modulus squared of the ratio is thus 
finite in any gauge theory. For a conformally-invariant gauge theory in d = 4, one finds 
the very simple expression 



r(Q2) 


2 r 


r(-Q2) 


= exp 



TT 



IK (a.(Q')) 



(2.16) 



Since all quantities in Eq. (2.16) have a precise nonperturbative definition, and since it 
provides a finite, unambiguous resummation of perturbation theory, Eq. (2.16) can be 

argued to be an exact result. It is easy to see that it agrees with the result for the four- 
loop form factor ratio for quarks in QCD [54], after setting all j3 function coefficients to 
zero. 
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3. Jet and Soft Functions in the Factorized Amplitude 



The general arguments for the factorization of QCD ampHtudes and cross sections, isolating 
the contributions responsible for long-distance singular behavior, have been known for some 
time [11]. In this section, we review the operator interpretation of this factorization, as 
described originally by Collins [47]. We then discuss the one-loop corrections to each of 
the relevant factors in dimensional regularization. 

At all orders, three basic physical principles apply to hard-scattering amplitudes, in- 
cluding the form factors under consideration. First, soft gluons decouple from hard virtual 
partons at leading power in the hard scale, since each soft gluon insertion in a hard subdi- 
agram adds to the diagram a new propagator far off the mass shell. Second, virtual hard 
partons collincar to an external hard parton effectively decouple from the remainder of the 
hard subdiagram, becoming insensitive to the energy and spin of fast partons moving in 
different directions. Finally, soft gluons decouple from jets, because their long wavelength 
does not allow them to discriminate features of a narrow jet other than its overall color 
and direction. 

To express an amplitude in factorized form unambiguously, we identify operator ex- 
pressions generating the leading contributions in each relevant region in momentum space, 
and make subtractions appropriate to avoid double counting. In the case of the form factor, 
the resulting factorization [47] is depicted in Figure 1 and can be expressed as 



n 



(3.1) 



For definiteness, in Eq. (3.1) we have in mind the time-like form factor for a massless quark, 
so that Pi=P2 = and {pi -\-p2)'^ = Q^- We also define the quark and antiquark velocities 
(5i via p^ = ((5/\/2)/3f ) so that (ii • j32 = 1, while the vectors define the directions of 
auxiliary gauge links to be discussed below. The hard function C summarizes the short- 
distance contributions to the form factor, and is finite as e — >^ 0. In order to define the 
remaining functions appearing in Eq. (3.1), it is useful to introduce first a notation for the 
Wilson line which describes the cikonal couplings arising in the soft and collinear limits. 
These couplings are generated by the operator 



^n(A2,Ai) = Pexp 



■A2 

'\g j dXn ■ A{Xn) 



(3.2) 



describing a gauge link in direction n'*. 

The 'partonic jets' J appearing in Eq. (3.1) are matrix elements for the transition 
from partonic states to the vacuum, mediated by the corresponding partonic (in this case, 
quark) field. Their operator definition is 

J (^^^,as{n^),e) u{p) = (0 |$„(oo, 0) V(0) \p) . (3.3) 
\ n jj, J 
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Figure 1: A graphical representation of the factorization of the form factor 

The spinor u{p) has been inserted to normahze the jet to unity at zeroth order; color and 
Dirac indices are implicit. The gauge link in the direction has a double role: on the 
one hand, it makes the matrix element gauge invariant; on the other hand, it mimics the 
coupling of gluons collinear to the incoming parton (say, a quark) to the opposite moving 
hard parton (say, the antiquark). Notice that, so long as 7^ 0, the function J (and 
similarly J, defined below) is invariant under rescalings of the vectors rij. Factorization 
makes it convenient to consider the vectors in Eq. (3.1) to be space-like [47]; a typical 
choice is n-i = Pi — P2 = —n2', the choice of vectors is otherwise free. This freedom can 
be used to derive the evolution equation (2.2) [6]. Clearly, J has infrared divergences, as 
well as collinear divergences associated with gluons collinear to p^. 

Next, we introduce the soft function S, as the vacuum expectation value of two light- 
like Wilson lines in the directions /3i and /32, 

5 (/?i • /?2, as(/^'), e) = (01^^2(00, 0)$;3i (0,-00) |0) . (3.4) 

The soft function S is the eikonal limit of the full form factor, and thus contains double 
poles at every order, associated with gluons that are soft and collinear to either hard leg. It 
has several useful properties. First, it is a pure counterterm in any minimal regularization 
scheme, because all its Feynman diagrams have no mass scale. In addition, because it 
is defined purely in terms of Wilson lines, it exponentiates according to the general non- 
abelian exponentiation theorem [55, 56], and its logarithm can be expressed in terms of a 
specific subset of Feynman diagrams ('webs'), before the soft loop momentum is integrated 
over. An important feature of the defining matrix element (3.4) is that it is invariant under 
boosts along the P1P2 axis in any frame where these velocities are back to back. On the 
other hand, while one might expect S, from its operator definition, to be invariant under 
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rescalings of the velocities, for light-like (3i this invariance is broken (as we will further 
discuss below); thus, S can depend on the velocities, through the combination /3i • (32 only. 

There are a number of subtleties in the evaluation of purely eikonal functions, including 
S. As a sum of scaleless integrals, S vanishes before renormalization on a diagram-by- 
diagram basis. This is perhaps why classic works on the renormalization of Wilson loops 
determine anomalous dimensions by introducing explicit mass scales [32, 33, 34], such as 
time-like lengths f3f > and/or cutoffs on the lengths of the Wilson lines $3-. This 
method was also used in Ref. [14]. Since here we are discussing specifically dimensionally- 
regularized amplitudes, we prefer to employ dimensional regularization throughout. To do 
so consistently, we must identify at each order the infrared-regularized coefficient of an 
ultraviolet pole. This is possible precisely because the logarithm of S has only a single 
ultraviolet divergence, which can be isolated systematically. 

As a practical matter, to evaluate dimensionally-regulated integrals it is convenient 
to rescale the light-like velocities /3j to have units of mass. In this way, standard shifts 
of loop momenta can be carried out. This can be done without loss of generality in the 
Wilson lines ^/j. that define S, by simple changes of variables. Below, we shall replace (3i 
by (jLt/\/2)/3j, with /3i • /32 = 1. Then, for the choice 11 = Q, the velocities are identified with 
their corresponding momenta. For now, however, we keep the scales of the /Jj arbitrary, 
and will continue to refer to them as 'velocities'. 

We now turn to the exponentiation of poles in S. CoUinear contributions to the soft 
function S can be factorized from purely soft contributions in the same manner as for 
the full form factor, but with partonic jets replaced by eikonal jets, for which we will 
provide definitions shortly. As a consequence, S satisfies an evolution equation analogous 
to Eq. (2.2), and the resulting solution is of the same form as Eq. (2.8) for the partonic 
form factor, but with — replaced by everywhere, and with an e-in dependent single- log 
function Gcik replacing G. The functions 7_ft"(as) and -fC(e, as), on the other hand, are the 
same as in Eq. (2.8), because the eikonal form factor matches exactly the partonic one in 
all infrared-singular regions, including the infrared-collinear ones. The arguments leading 
to Eq. (2.11) still hold, and we are led to 

2j^~^ '^"'^ (^1 • /?2, a (^^ e) j 
-\lK{a (^^e)) log 



(3.5) 



Notice that the lack of explicit e dependence of Geik ensures that iS is a pure counterterm, 
consistent with its diagrammatic interpretation. We expect the real part of the function 
Geik, computed in back-to-back kinematics, to be related to the anomalous dimension Fdy, 
defined and computed at two loops in Ref. [50]^. 



The "leading tramscendeiitality" term in r^Y: proportional to C(^)) controls the single poles of the 
polygonal Wilson loop expectation values computed recently at two loops in = 4 super- Yang-Mills 
theory [45]. 
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Including in the factorization both partonic jets and the full eikonal form factor clearly 
double counts the soft-collinear regions. This can be avoided if one divides by eikonal 
versions of the two jets, which are defined as 

= (0|$„,(oo,0)$;3,(0,-oo)|0) , 

= (0|$;3,(oo,0)$„,(0,-oo)|0) . (3.6) 

These jets are also pure counterterms in dimensional regularization, because they do not 
depend on any mass scale. They have soft-collinear enhancements from gluons moving in 
the /3j directions, matching those of both the soft function S or the partonic jets J. As a 
consequence the ratio J/ J' has only single collinear poles at every order, associated with 
hard collinear radiation, while the ratio 

^ 2n x _ 5(/3i •/32,as(M^),e) 

has only single infrared poles associated with soft gluons emitted at wide angles from the 
hard partons. In Eq. (3.7) we have noted that the function iS can only depend on the 
homogeneous ratio 

p... (-ft-fe)%;.^ (3.8) 
{-Pi ■ ni) {-(32 ■ n2) 

The reason is that invariance under separate rescalings of the velocities /3i, which was 
broken for 5, must be recovered in <S, which contains the complete (3 dependence of the 
form factor; furthermore, homogeneity in rii is built into the eikonal Feynman rules. The 
simplicity of these relations is a direct result of eikonal exponentiation. In the exponent, 
collinear and soft regions enter additively and universally [57]. 

The one-loop diagrams associated with the functions entering Eq. (3.1) are easily 
evaluated using eikonal Feynman rules where appropriate. We give the results below in 
the MS scheme; thus all factors of log(47r) and 7e are absent, having been absorbed into 
the definition of the renormalization scale ^. 

The soft function receives a one-loop contribution only from the vertex correction 
diagram, since self-energies on eikonal light-like lines vanish like (3f. The full one-loop soft 
function is then given by the UV counterterm for the time-like vertex correction [58]^, 



<5« (/3i-/32, e) = -^CF 
47r 



^--log(-/3i-/32) 



(3.9) 



^Note that the coefficient of — /3i • 02 in the argument of the logarithm in Eq. (3.9) can be changed 

by rescaling the eikonal Feynman rules. Associating with each gluon emission, for example, a factor of 
K;/3^/(fv/3 • fc), instead of the usual factor /3^/(/3 • k), rescales the argument of the logarithm by a factor 
K^. This ambiguity is associated with the broken invariance of the function >S under rescalings of and 
corresponds to a choice of scheme in the renormalization of <S, which was discussed above. Once again, 
this ambiguity does not affect physical quantities: the dependence on k cancels between the soft function 
>S and the eikonal jets J , as discussed below Eq. (3.8). Notice also that the invariance under rescalings of 
the vectors nf in the jet functions is not broken, since ^ 0, so that there are no collinear divergences 
associated with them. 
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Matching Eq. (3.5) and Eq. (3.9), we recover, as expected, that 7^^'' = 2Cf at one loop, 
while^. 

(A ■f32) = --CF log (-A • 132) . (3.10) 
vr 

The eikonal jet J receives contributions at one loop from both the eikonal vertex correction 
and the self-energy diagram on the eikonal line along the direction. This eikonal self 
energy, J^2 below, is a single pole pure counterterm at one loop, which is common to the 
partonic jet J and to the eikonal jet JT". It cancels in their ratio, but contributes to their 
respective evolution equations. The vertex correction to J , which we denote by J\fi also 
contributes only through a counterterm. The complete one-loop eikonal jet is then 

7(1) - 1 7(1) . 7(1) 

1 1, / 



r{l) _ as ^ I 



e2 e ^\2{-(3-n) 



= -Cf-. (3.11) 

^ 27r e ^ ' 

As expected, the soft-coUinear double pole of the eikonal jet is one half of the corresponding 
pole in the eikonal form factor S. Here and below, the one-loop self-energy counterterm 
(for the n eikonal) is multiplied by 1/2, which reflects the removal of the square root of the 
residue of the relevant two-point function in a normalized S-matrix element. 

Turning to the partonic jet J, we encounter in its vertex correction a one-loop diagram 
that is not simply a pure counterterm. The full one-loop result is the sum 

^^^^ = ^^^^^ + ^^^^ + ^4^ (3.12) 
with J^^ given as above, while Jy ^ is the quark-eikonal vertex correction, given by 



.(1) ( jP ■ \ _ as 



, , 2 + log . 

+ 1 + '"^ ((4^) + ^ + + o (3-") 

finally, Jp^ is the pure counterterm self-energy on the quark leg, 

4'He) = ^CF'. (3.14) 

Note that, as expected, the vertex correction Jy ^ gives the same double pole as the eikonal 
jet, Eq. (3.11). 



■* Equation (3.10) is consistent with Refs. [50, 59], which find that F^^l = F^l^ = 0, provided one chooses 

the subtraction scheme for coUinear poles corresponding to k = 1 of the footnote above. In this scheme the 
argument of the logarithm is — /3i -02, as shown in (3.9), so that G^^^ is purely imaginary for back-to-back 
time-like kinematics, giving a vanishing contribution to one-loop cross sections, while it vanishes for space- 
like kinematics. Note also that the "DIS" contour used in Ref. [59] differs from the space-like configuration 
considered here, so CiH has no correspondence with the defined there. 
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Collecting all the ingredients, and applying Eq. (3.1), we expect to reproduce all the 
infrared and coUinear poles of the form factor at one loop. Indeed we find 



r(i) - c{i) (a R 7(1) ( ■ "1) 

pole \^^' -^"{Pi- P2, e) + Jy p^ig 



J. 



(1) 

V, pole 



{P2 ■ 112) 



V 



a 



47r 



2 2 / 



V 



6) +4^) (e) 



(3.15) 



which matches the pole structure of the one-loop form factor in dimensional regularization, 



r(i) 



Equation (3.16) also implies that 



r(l-2e) 



TT 



Cf 



log 



(3.16) 



(3.17) 



in agreement with Ref. [7]. Our task is now to construct an all-order expression for the 
function G in terms of the anomalous dimensions of the various functions building up the 
form factor according to Eq. (3.1). 

4. From Factorization to an Operator Intepretation for G{as) 

Let us begin our investigation of G{as) by considering the rcnormalization properties of 
the various functions entering Eq. (3.1). The partonic jet function J and the short-distance 
function C are multiplicatively renormalizable, with anomalous dimensions depending on 
the coupling ag but not on the infrared regulator e. Eikonal functions such as S and 
on the other hand, require extra care: in general, their anomalous dimensions need 
infrared regularization in their own right, because of the overlap of collincar and ultraviolet 
divergences in the renormalization of any 'cusp' singularity involving light-like Wilson 
lines [34]. 

Consider first the soft eikonal function S. The observation that it is a pure counterterm, 
of the general form of Eq. (3.5), leads to 



/^^log>S (/3i •/?2,Q;s(/x^),e) = ^(e, a^) — log<S (/?i • /32, ^^(m^), e) 

= - IS • /32,Q!s(/^^),e) , 



(4.1) 



where the singular anomalous dimension 75 is related to the function Geik and to the cusp 
anomalous dimension 7^ by 

75 {Pi ■ P2, as{ii^),e) = - Geik {Pi ■ P2, ocgin^)) + 2 Jq i^i^^^ ^)) 

= - Geik • P2, oisili^)) - K (e, as^li")) . (4.2) 
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Here, in the second equality, we have used Eq. (2.10). We observe that 75 has a single 
infrared pole determined by the cusp anomalous dimension ^k- In order to work with 
infrared-finite anomalous dimensions, one can instead consider the function <S, defined in 
Eq. (3.7), which carries only single infrared poles due to wide angle soft radiation. It obeys 



H^^ogS {pi2,ots{li^),e) = - 75 (pi2,as(/x2)) 



(4.3) 



Finally, one may define 



^ 1 7 
/X— log J 

djjL 



{Pi ■ njf 



djjL 



-lj{0is) , 

-7c(pi2,as 



(4.4) 



where the functional dependence of jc on /x is dictated by the requirement that the form 
factor as a whole not be renormalized, which implies, using Eqs. (3.1) and (3.7), 



Is (Pi2> "s) + 7C {Pi2,as) + 27J (as) = . 



(4.5) 



At one loop, we can compute 7j by combining the terms specified in Eq. (3.12); similarly, 
75 is derived from the one-loop results for S and J, combined as in Eq. (3.7). We find 



(1) /X «s ^ 

7-^ (P12) = —Cf 



7?' 



3 

TT 4 



Cf , 



(4.6) 



where we note that 7^^ equals the one-loop anomalous dimension of the quark field. We 



then derive for 7^^ using Eq. (4.5), 

7c^ (P12) = 



TT 



Cf 



Mr 



(4.7) 



Having exhibited these one-loop examples, we continue with the general discussion. 

The next step is to consider the dependence on the vectors nf, which enter the form 
factor through the jet functions J and their eikonal counterparts J'. Following the reasoning 
of Rcfs. [4, 5] in axial gauge, generalized to arbitrary gauges in Rcf. [6], we begin by 
observing that the form factor must be independent of nf. Defining Xi = {—Pi • rii)^ /n?, 
we can write 



OXi 



,as(//^),e = 



(4.8) 



When we apply this consistency condition to the factorizcd cross section in Eq. (3.1), 
derivatives with respect to of the partonic jet functions Jj separate into ultraviolet- and 
infrared-dominated terms, according to 



d d d 

Xi ^ log Ji = -Xi — logC + Xi — log Ji 

_ 1 

^ 2 



Gi[xi,as{iJ. ),e) +/C(as(/x ),e 



(4.9) 
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where the second line defines the sum of the functions Qi and /C. As is clear from this 
definition, and in analogy with Eq. (2.2), the function Qi carries the scale dependence, 
but is finite as e — > 0, while /C is a pure counterterm. At one loop we find, directly from 
Eq. (3.12), 

^f^(..e)^gc.(log^ + l)+0(.), 

/CW(e) = ^Cf^-. (4.10) 

We can now relate the functions Qi and K. to matrix elements of fields in the presence of 
Wilson lines. In fact, both partonic and eikonal jets depend on and on the velocity 
only through the combination x = {—(3 ■ nf' /n^. One can thus simply relate their x 
dependence, given in Eq. (4.9), to their dependence, using 

dJ V? dJ 

p-n- = 4.11 

op ■ n p-n O'nr 

and similarly for J, with p'* replaced by (3^. Prom the definitions of the jet functions, 
Eqs. (3.3) and (3.6), and from the behavior of an ordered exponential under variation with 
respect to the curve, we readily find expressions for the derivatives of the jet functions in 
Eq. (4.9), extending results found by Collins in QED [47]. 

For both partonic and eikonal jets, a derivative with respect to the vector replaces 
an ordered exponential in the direction by the integral of a field strength over the 
original path, sandwiched between Wilson lines. Written explicitly, we can use this result 
to determine the p ■ n dependence of the jet functions (3.3) in terms of matrix elements, as 



d 

p ■ n 



dp ■ n 



( 7) ' TZ^^ 
log J ( 2-^'««(/^^^ 



.-1 t/ ^ /n, 9 



J- 



p • n 2pQ 



uHp) (0|— $n(oo,0)V(0) \p) 



p ■ n 



1 

J~' TT- u\p) / d\ A (0| $„(oo, A) p^ 7f (ig F^,{\n)) $„(A, 0) V(0) \p) 
^Po Jo 

= ^(^^2^'«^(/^')'e) +^(«s(/x'),e) , (4.12) 

where the final equality is simply a restatement of Eq. (4.9). Similarly, as the first equality 
in (4.9) makes clear, /C is found directly from J' as 

„2 /-oo 

= - — J-W dA A (0|$„(oo, A) ^^rf (i^ F^.(An)) $„(A, 0) $^3(0, -00) |0) 
P ''f^ Jo 

= }C{as{n^),e) . (4.13) 

In the perturbative expansions of these matrix elements, the field strength operator pre- 
vents an unphysical gluon from coupling to the eikonal. At each order, the corre- 
sponding vertex cannot appear in any jet-like subdiagram that provides a coUinear pole. 



- 16 - 



It can, however, appear in subdiagrams that carry ultraviolet and infrared momenta [4], 
and, in the latter case, are associated with an infrared pole in dimensional regularization. 
The inverse jet and eikonal-jet factors multiplying the matrix elements cancel the residual 
coUinear singularities, which factorize. The remaining terms give, order-by-order, the Q 
and /C functions above, from the short- and long-distance non-coUinear regions to which 
the field strength vertex can contribute. The function Q is seen to be the difference of 
two gauge-invariant matrix elements, both involving the field strength and Wilson lines, 
derived from, and normalized by, the partonic and eikonal jet functions. 

It is worth noting that the — limits of Eqs. (4.12) and (4.13) are singular. In this 
limit, the matrix elements become boost-invariant, and at the same time develop coUinear 
singularities for gluons in the direction. As for the soft function S, the jet functions, 
which are sums of pole contributions only, are no longer scale invariant in the light-like 
vectors f3 and n. 

We can now turn to the determination of the Q dependence of the full form factor, 
Eq. (2.2). r depends on Q directly through the short-distance function C, and indirectly 
through the partonic jets Jj, which depend on external momenta through Pi-rii. Prom the 
factorized expression, Eq. (3.1), using Eq. (4.5), one easily derives 



Because the pole terms /C(e, a^) are independent of the kinematic variables, they are equal, 
and we have K = 2IC. Comparing Eq. (4.14) with the original evolution equation, Eq. (2.2), 
we finally find an expression for G, in terms of the anomalous dimensions of the soft and 
eikonal jet functions, and in terms of the functions Qi, defined by the matrix elements of 
Eqs. (4.12) and (4.13). Explicitly, we have 



which is easily verified to hold at one loop, making use of the results of Section 3. 

5. Eikonal and CoUinear Contributions to Subleading Poles 

The single poles of form factors are generated by functions G{as), one for each type of 
parton, which have been related to gauge theory matrix elements via Eq. (4.15). In this 
section, we will link G{as) to two anomalous dimensions, one stemming from the coUinear 
evolution of parton distributions, the other from the eikonal form factor, S. A relation 
following from this structure was verified empirically to three loops in Ref. [49], based 
on an earlier observation of Ref. [48] . A similar connection has been established at finite 
order between single-logarithmic contributions to the Drell-Yan cross section and coUinear 
evolution kernels in Ref. [17]. These empirical observations are established here to all 
orders in perturbation theory, exploiting a connection between the form factor and parton- 
in-parton distributions, which follows from factorization and which was noted already in 




(4.14) 




(4.15) 
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Ref. [60]. To derive our result, we first have a look at the analog of the form factor in 
parton evolution. 



5.1 Factorization for Virtual Contributions to Pjirton Distributions 

Let us begin by considering the standard definition of the light-cone distribution for a 
parton of flavor i, carrying momentum fraction x, in a parent parton of the same flavor. 
For a quark, for example, one writes 

^) = 4^ / ^e-^^-^'-^(p|V^,(A/3) 7 • /? MX, 0) V,(0)|p) , (5.1) 

where p is the momentum of the parent quark, which we can take, say, along the (+) 
direction; /3 is then an auxiliary light-cone vector along the (— ) direction, and the Dirac 
projector 7 • /? selects the relevant components of the quark field, while the Wilson line 

ensures gauge invariance. An analogous definition applies for the gluon-in-gluon distri- 
bution, with the Wilson line in the adjoint representation. Note that there is no explicit 
p ■ (3 dependence in the parton distributions when they are defined in this boost-invariant 
fashion in the MS prescription. 

In order to single out the virtual contributions to the parton distribution, we proceed 
as follows. The gauge link can be split by extending it to light-like infinity along the /3 
direction, according to 

$^(A,0) = $;3(A,oo)$;3(oo,0); (5.2) 

one can now insert a complete set of states between the two Wilson lines, and then iso- 
late the contribution of the vacuum. This gives the virtual contribution to the parton 
distribution at the amplitude level as the correlator 

r,/, [P—^,as{li^), = (0|$;3(oo, 0) V,(0) \p) , (5.3) 

coupling a single-particle state to the vacuum through the action of the partonic field ijjq 
and of a gauge link in the same color representation. We shall define this matrix element 
as a sum of pure pole terms, consistent with its interpretation as part of an MS parton 
distribution function. 

Clearly, the amplitude ^q/q in Eq. (5.3) is closely related to the partonic jet, J, in 
Eq. (3.3). In fact, the only difference is that the gauge link is now in a light-like direction 
opposite to the parton momentum. As a consequence, F^/g can be factorized in the same 
manner as the full partonic amplitude, Eq. (3.1), into short-distance, jet, and soft functions. 
Now, however, we need a separate partonic jet only for the incoming line, since the coUinear 
singularities of the outgoing gauge link match the collinear singularities of the soft function 
S in Eq. (3.1). We can then write 



(5.4) 
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where we have introduced the velocity four-vector /3p associated with the momentum p. 
The function Cj is a short-distance coefficient chosen to cancel all terms that are finite for 
e ^ in J, because Tg/g is defined as a sum of pole terms only. We are assured that the 
function Cj exists, because of the exponentiation of all pole terms. Using the fact that S 
and are also pure pole terms, we may write. 



2-^,Ols{fi ),€ 



(5.5) 



pole 



Because ^q/q is n-independent, all n-dependence in poles on the right-hand side of Eq. (5.4) 
is guaranteed to cancel in the ratio of jet functions, leaving only finite n-dependence, which 
is cancelled by Cj. This cancellation is possible simply because coUinear singularities are 
independent of n'' [4], leaving only soft contributions, whose (exponentiating) poles match 
between J and J', and short-distance contributions, which are cancelled by Cj. We observe 
that the ratio of partonic to eikonal jets in Eq. (5.4) is the same as in the basic factorized 
form, Eq. (3.1). 

Inserting Eq. (5.5) into Eq. (5.4), Tg/g can be represented as. 



pole 



(5.6) 



where we have chosen the magnitude of the four- velocity f5p so that p = (/i/\/2)/3p, and the 
scalar product p- (3 = {n/\/2) (3p ■ (3. Thus /i is the only remaining scale, which appears only 
as the argument of the coupling. As usual, Eq. (5.6) involves exponentiating double poles 
from S, which cancel when combined with the real emission contributions to the parton 
distribution, leaving behind only single, collinear poles that define the splitting functions. 

Now, following Ref. [60], we introduce the eikonal counterpart of the correlator (5.3), 
which will represent the virtual contribution to an eikonal parton distribution. This is 
naturally defined as the soft function S itself. 



(5.7) 



which, once again, is an exponential consisting entirely of pole terms. We can now ex- 
ploit the fact that MS parton distributions can be defined, in moment space, simply as 
exponentials of the integrated collinear anomalous dimension, as was done in Ref. [60]. 
Furthermore, the eikonal approximation is accurate for real final-state radiation, up to in- 
verse powers of the Mellin variable N. The ratio of the virtual contribution for the parton 
distribution to its eikonal counterpart must thus be given entirely by the virtual term of 
the corresponding diagonal splitting function, B^^\as), whose normalization is defined by 



Pii{x) 



1 



x) + o((l-x)°) . (5.8) 
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Taking the ratio of Eq. (5.7) and Eq. (5.6) we then find, 



J 



pole 



cik 



J ((^, asin^),e) T^/, iPp ■ P, a,(/x2), e) 



exp 



(5.9) 



where the second relation^ mirrors the results obtained in Ref. [60] for the complete par- 
ton distribution. At one loop, B^^\as) = (3/4) Cf (os/tt) for quarks, and B^^\as) = 
{bo/4:) {ag/ir) for gluons. Prom Eq. (5.7), the eikonal vertex T^^'^ and S are to be computed 
in the same fashion in perturbation theory, identifying their eikonal velocities with mo- 
menta (/x/-v/2)/3p and (/i/\/2)/5, for incoming and outgoing lines, respectively. Consistency 
in the factorization formula Eq. (3.1) then requires the same treatment of the velocity Pp 
in the function J on the left-hand side of the first equality in Eq. (5.9). 

The result in Eq. (5.9) expresses the purely collinear single poles of the ratio between 
the partonic jet function and its eikonal counterpart in terms of the virtual contribution 
Bs{as) to the splitting kernel of the appropriate parton flavor. We now use this result in 
conjunction with our basic factorization formula, Eq. (3.1), to get another simple expression 
for the function G. 



5.2 Relating Form Factors to Collineeir Evolution Kernels 

We proceed by using Eq. (5.5), followed by Eq. (5.9), in Eq. (3.1) for each jet, and exploiting 
renormalization-group invariance of the full form factor to set /x^ = Q^. This expresses 
r((5^,e) in terms of the eikonal soft function S, the virtual evolution kernel Bs, and the 
finite factors in the partonic jet functions, as 

T{l,as{Q%e) = c(l,^-^^,as{Q%e) iflCj (^^^ ,a,{t^%e)V 



X exp 



10 



Bs{a{e,e)) 



= C{as{Q^),e) 

de 

10 e 



X exp 



2 Bs{a{e,e)) 



S{Pi-P2,as{Q^),e) , 



S{Pi-P2,as{Q^),e) , 



(5.10) 



where in the second relation we define the function C{as, e) to include all factors that are 
finite at vanishing e. In this function, all dependence on jet directions cancels. 

In Eq. (3.5), wc have an exponentiated form for the soft function S, but to make direct 
contact with the standard form factor notation, Eq. (2.11), in terms of G and 'Jk, we 



^ There is a factor of 1/2 in the exponential of Eq. (5.9) relative to the corresponding relation in Ref. 
[60], because we are computing an amplitude here rather than a parton distribution. 



- 20 - 



need an expression for C as well. This expression can be simply obtained, using again the 
integrability of the coupling for e < 0. We write 



C 



,as{Q'^),€ ] = exp 



= exp 



e 

1 /■^'d?^ 

2 Jo C 



2 (^^^) 



, a 



(5.11) 



In the definition of Gq we have used the independence of the form factor on the choice of 
the eikonal vectors rii, and have inserted a factor of 1/2 to conform with the normalization 
of Eq. (2.8). We note that although G-^ is finite at e = 0, the integral in Eq. (5.11) 

produces no poles, since the logarithmic derivative with respect to the scale generates a 
positive power of e. Inserting Eq. (5.11) and Eq. (3.5) in Eq. (5.10), we find our final 
exponentiated result. 



,as(M^),e = exp 




Geik(l,a(e^e))+2B^(a(e^e)) 



+ G^ {a{e, e),e)-^-jK («(C', e)) log } , (5.12) 

where we have set /3i • /32 = 1 in Geik- Comparing this result with Eq. (2.11), we conclude 
that 



G{l,as,e) = 2Bs (as) + G^ik (l,"^) + G^{as,e) , 



(5.13) 



The function G is thus the sum of three terms: twice the coefficient of (5(1— x) in the relevant 
parton splitting function, the single-logarithmic anomalous dimension of the eikonal form 
factor, and finally a term associated with the running of the coupling in the infrared- 
finite hard-scattering function. The latter term is proportional to the d-dimensional beta 
function, Eq. (2.4) (see Eq. (5.11) and the Appendix); hence it vanishes as e — in a 
scale-invariant theory. 

Finally, comparing Eqs. (4.15) and (5.13) for G{as), and referring to Eq. (4.12) and 
Eq. (4.13), which relate Q to nonlocal matrix elements involving the field strength, we 
find that the moment-independent term in the evolution kernel is determined by the same 
matrix elements of the field strength, and by a combination of anomalous dimensions 
of eikonal and local operators, including the logarithmic derivative of the functions Cj 
introduced in Eq. (5.4). Explicitly, 



Bs 



=1 



/3(a,)5^1ogC.. 



1 



Geik 



1 



^75 - 7J ■ 



(5.14) 



Although individual terms in this new relation for the function Bs{as) depend on rii and 

Pi, this dependence cancels in the sum. 

Once again, Eq. (5.13) can easily be tested at one loop using the results of Section 3. It 
can be tested further, up to the three-loop level, by comparing with the results of Ref. [49], 
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where the relation between the function G and the virtual splitting kernel was emphasized. 
Indeed, Eq. (20) of Ref. [49] shows that up to three loops the perturbative coefficients G^''^ 
of the function G are given by the sum of 'maximally non-abelian' terms f^''\ plus twice the 
virtual splitting kernel contributions Bg^\ as in Eq. (5.13), plus remainders proportional 
either to e or to the (3 function coefficients. In the Appendix we show that such terms are 
precisely the ones that arise in the expansion of a function of the coupling that is defined as a 
total derivative with respect to the scale, as is the case for G^j. It is natural then to identify 
f(^) with the corresponding perturbative coefficient of G^ik- Indeed, f^^^ = 0, consistently 
with Eq. (3.10) in the MS scheme and for space-like kinematics. Furthermore, one easily 
verifies that, when brought to the same normalization, /^^^ is one half of T^l = T^y^ as 
computed in Ref. [50]. The factor of 1/2 is expected, since /^^^ contributes to an amplitude 

(2) 

while contributes to a cross section. 

Equations similar to Eq. (5.13) have appeared in the description of the anomalous 
dimensions of effective currents (or their matching coefficients) in soft collinear effec- 
tive theory. In Ref. [18] it was noticed that the empirical relations found through three 
loops [17, 48, 49] imply that the subleading-logarithmic part of the effective-current anoma- 
lous dimension, denoted there Bi{as), is given just in terms of Bs{as) and the function / of 
Ref. [49] ; the terms drop out of the anomalous dimension. In Ref. [20] a similar relation 
was found, and / (called 7^ in Ref. [20]) was identified with the anomalous dimension for a 
momentum-space Wilson loop associated with the Drell-Yan process. Analogous relations 
were obtained for the Drell-Yan process in Ref. [61] and for deep-inelastic scattering in 
Refs. [19, 62]. 

The identification of the coefficients f^^^ with the eikonal (Wilson-line) quantities G^^ 
neatly explains a couple of their properties found empirically through three loops [48, 49] : 
the relation fg'^^ /Ga = fq'^ /Cp, and the maximally non-abelian color structure of these 
quantities. The non-abelian exponentiation theorem for eikonal graphs [55, 56] implies 
that only "color-connected" graphs composed of single gluon webs (along with fermion- 
loop insertions) contribute to Geik- Through /c = 3 loops, all such graphs have color factors 
of the "maximally non-abelian" form CjC^^'^^rij, where Cj is the Casimir factor for the 
eikonal line, Cp for quarks and Ga for gluons. This form breaks down at four loops, due 
to the existence of color factors that cannot be expressed in terms of quadratic Casimir 
operators, as in the case of the four-loop beta function in QCD [63]. 

As mentioned above, in a conformal theory there is no contribution to Eq. (5.13) from 
G^ as e — 0, so we have 

G(l, a,, 0) = Geik (1, as) + 2Bs (as) . (5.15) 

The eikonal quantity Geik carries no information about the spin of the parton, only its 
color (representation under the gauge group). Thus the spin-dependence of G(1,q;s,0) is 
all carried by the virtual part of the splitting kernel, Bs{as). Many conformal theories are 
supersymmetric: in this case, if two partonic states belong to the same supersymmetry 
multiplet, then they are in the same gauge-group representation, and Geik is the same for 
both. The leading-twist operators whose anomalous dimensions yield Bs{as) will also be 
related by supersymmetry. By Eq. (5.15), the values of G(l,as,0) for these states should 
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be the same too. This result can also be seen via supersymmetry Ward identities which 
relate the appropriate S'-matrix elements [64], and thereby imply that the corresponding 
single 1/e poles have to be identical. It would be interesting to see if the simple compound 
representation of G(l, as, 0) in Eq. (5.15) can help in the study of its properties in A/" = 4 
super- Yang-Mills theory. 

6. Concluding Remarks 

We have reviewed the resummation of poles in dimensionally-regularized singlet form fac- 
tors for QCD and related theories, using the factorization properties of their amplitudes. 
Revisiting the basic evolution equation and its solution, we observed the strikingly simple 
connection between the analytic continuation of the form factor in a conformal theory and 
the cusp anomalous dimension, given by Eq. (2.16). 

In view of of the universal nature of coUinear poles in dimensionally-regularized ampli- 
tudes, Eq. (1.1), we have investigated the origin of subleading poles in the form factor. Our 
analysis extends the familiar relationship between the leading poles and the cusp anomalous 
dimension. We have determined the origin of the non-singular function G{as) in the basic 
evolution equation, Eq. (2.2). Relying on the operator content of the factorized jet and 
soft functions, we found Eq. (4.15), which relates G{as) to the short-distance function C, 
and to 7^ and jj, the anomalous dimensions of soft and jet functions respectively, as well 
as to matrix elements involving the field strength and Wilson lines. The structure of this 
result is made more transparent in Eq. (5.13), which expresses G{as) in terms of two uni- 
versal quantities (the single-pole anomalous dimension of the eikonal form factor, and the 
coefficient of (5(1 — z) in the diagonal partonic evolution kernel), plus a process-dependent 
short-distance contribution, which generates no singularities and is proportional to e in 
the conformal limit. The same pattern has been noted in deep-inelastic scattering and the 
Drell-Yan process [17, 18, 19, 20], and, through the relation between form factors and am- 
plitudes, it will appear in subleading logarithmic corrections to any threshold-resummed 
jet cross section [65]. Equating our two expressions for G{as) provides an interesting new 
relation between Bs{as), other anomalous dimensions, and matrix elements of the field 
strength. 

In summary, unlike leading poles, nonleading poles in form factors and fixed- angle 
scattering amplitudes have a compound structure, even in conformal theories. Nevertheless, 
all contributions to G{as) that generate infrared poles have a well-defined and universal 
origin, in terms of matrix elements and anomalous dimensions in the massless gauge field 
theory. In the context of planar J\f = 4 super- Yang-Mills theory, the explicit operator 
interpretations for G{as) may aid efi^orts to use integrability to determine its exact coupling- 
constant dependence. 
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A. Compcirison to Fixed-Order 

The form factors for quarks and gluons in QCD have been evaluated recently at three 
loops in Refs. [49, 54], allowing for a stringent test of our result, Eq. (5.13). Ref. [49] 
(MVV below), in particular, gives an explicit expression for the single-logarithmic function 
G{as, e) up to three loops, as a sum essentially in the form given in Eq. (5.13) above. In 
this sum, the function that we denote by G^(as, e) corresponds to a set of terms containing 
the coefficients in Eq. (20) of MW. Similarly, Geik can be identified with terms labeled 
/ in MVV, whose universality is noted, without an explicit discussion of their origin. 

The key feature of Eq. (20) in MVV is that the terms in the function G that are not 
accounted for by the virtual splitting kernel or by the eikonal function / are proportional 
either to e, or to the coefficients of the /? function. We want to verify that all these terms 
are precisely of the form that follows by requiring that they are coefficients in the expansion 
of a total derivative with respect to the scale of the running coupling, as is the case for our 
function G-^{as,e), Eq. (5.11). 

Let us begin by working out the consequences of our definition of G-^{as, e), which is 
of the form 

Gc{as{fi^),e)=2fi^^E{asiii^),e) , (A.l) 

for some function E{as{n'^),e) that is an expansion in both as(/x^) and e, and which is 
finite at e = 0. In terms of our factorization analysis, E{as,e) is simply the logarithm of 
the finite coefficient function G in the factorized form factor, in the scheme in which the 
coefficient function is defined to absorb all finite terms in the expansion of the partonic jet 
functions. Expanding G^ and E in powers of and e, we write 

oo oo oo 

n=l n=l m=0 

and similarly for E. Next, we use the fact that E depends on the scale only through the 
coupling as(^^). As a consequence, we can easily find an expression for the perturbative 
coefficients of Gq in terms of those of E. Using our normalizations for the /3 function, 
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Eqs. (2.4) and (2.5), we find 

2 



^ n— 1 

G^^\e) = -2neE(")(e) - - ^ S('=)(€) . (A.3) 



fc=i 

Since all ^('^^(e) are finite as e ^ 0, this implies 

n-l 

2 



^ n— 1 

GffiO) = "2 E ^&n-fc-i^^')(0). (A.4) 
fc=i 

We observe from Eq. (A.3) that the perturbative coefficients of Gq are sums of terms 
that are proportional cither to e or to the coefficients of the /3 function, as expected. 
Furthermore, combining Eq. (A.4) with Eq. (A.3), it is clear that one can determine G^Vo) 

recursively in terms of G^\e), with k < n. To match the notation of MVV, we proceed 
by defining 

g{e) = ^ [g{e) - g{0)] , (A.5) 

for any function g{e) with a finite limit as e — 0. The recursion starts with G^\o) = 0, 
so that Eq. (5.13) gives 

G(^) = 2S(^) + G^ + edff , (A.6) 

matching the one-loop Eq. (20) of MVV. 

Proceeding recursively, it is easy to see that at two loops one can write 

Gg)(0) = -|^a)(0) = ^G«(0). (A.7) 



?g^(e) = Gp(0) + eG'§\ 
Eq. (A.7) matches the two-loop result in Eq. (20) of MVV, 



Using G^ (e) = (0) + eG^ , and taking into account the different normalizations, 

O C-^ o 



G(^) = 25f + GS + |G«(0) + eGg'. (A.8) 
A short calculation yields also the three-loop expression 
Gg)(0) = -6o£;(^)(0)-|£;«(0) 

which again matches the /^-function terms in the three-loop result in Eq. (20) of MW, 
provided the normalizations are taken into account. 

An alternative way of expressing the solution of this recursion problem, without making 
use of the subtraction in Eq. (A.5), is to expand explicitly the coefficients G^\e) in powers 
of e, as done in Eq. (A. 2). It is straightforward to express the solution, at any order, in 
terms of the coefficients G^^ . One finds 

C , m 

4",:" = - leg;-" + ^Gg;-" + |Gg-> - "^G^f + . . . . (A.10) 
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Equation (A. 10) was derived in Ref. [53], as a solution to the problem of finding a function 
G{as{lJ'^)), independent of e, but capable, upon integration over the scale, of matching 
the simple poles of the form factor, generated by the function G{as{iJ.'^), e). The result of 
Ref. [53] can now be rephrased by stating that the function G(q;s(/x^)) must be given by 
G{as{n^),e = 0) plus the total derivative with respect to the scale of a finite function of e. 
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